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In this work we present an alternative way to look at electron diffraction in a transmission electron 
microscope. Instead of writing the scattering amplitude in Fourier space as a set of plane waves, we 
use the cylindrical Fourier transform to describe the scattering amplitude in a basis of orbital angular 
momentum (OAM) eigenstates. We show how working in this framework can be very convenient 
when investigating, e.g., rotation and screw-axis symmetries. For the latter we find selection rules 
on the OAM coefficients that unambiguously reveal the handedness of the screw-axis. Detecting the 
OAM coefficients of the scattering amplitude thus offers the possibility to detect the handedness 
of crystals without the need for dynamical simulations, the thickness of the sample nor the exact 
crystal structure. We propose an experimental setup to measure the OAM components where an 
image of the crystal is taken after inserting a spiral phase plate in the diffraction plane and perform 
multislice simulations on a-quartz to demonstrate how the method indeed reveals the chirality. 

The experimental feasibility of the technique is discussed together with its main advantages with 
respect to chirality determination of screw axes. The method shows how the use of a spiral phase 
plate can be extended from a simple phase imaging technique to a tool to measure the local OAM 
decomposition of an electron wave, widening the field of interest well beyond chiral space group 
determination. 

PACS numbers: 61.05.J-, 03.65.Vf, 41.85.-p 


I. INTRODUCTION 

An object is chiral when it is not superposable onto 
its mirror image by rotating and/or translating it. Two 
chiral objects that differ solely by mirror operations are 
called enantiomorphs and are labeled as the right- and 
the left-handed variant of the enantiomorphic pair. Chi¬ 
rality occurs from the macroscopic scale, e.g., a hand, 
to the nano-meter scale, e.g., molecules and crystals. In 
the latter it manifest itself not in the physical proper¬ 
ties of the material (e.g., density, melting point) but in 
the interaction with other chiral objects, such as circu¬ 
larly polarized light or other chiral molecules. This is 
of great importance in, e.g., the pharmaceutical industry 
since biological processes mostly are chiral. Enantiomor¬ 
phic compounds in drugs can have distinct effects on the 
human body, and it is essential for the drug to be enan- 
tiomorphically pure in order to have the desired effect. 
In crystals the chirality can be seen directly from the 
crystal’s space group lacking an improper symmetry el¬ 
ement (mirror, mirror glide plane, inversion center, or 
roto-inversion center). When on top of that the space 
group possesses a screw-axis apart from the 2i screw- 
axis, the space group itself is chiral as well. Namely, 
mirroring a crystal with a chiral space group, (see Table 
I), will change the space group of the crystal to the other 
one in the enantiomorphic pair. 

Determining the handedness of a crystal in a trans¬ 
mission electron microscope is a difficult task. In con¬ 
ventional transmission electron microscopy (TEM) the 
crystal is seen as a two-dimensional projection on the 
viewing screen and one can not distinguish between enan¬ 
tiomorphs from one single image, simply because mirror¬ 
ing the crystal along the projection plane changes the 


1. 

(P4i,P43) 

5. 

(P6i,P65) 

9. 

(P3i,P32) 

2. 

(P4i32,P4332) 

6. 

(P62,P64) 

10. 

(P3il2,P32l2) 

3. 

(P4i22,P4322) 

7. 

(P6i22,P6522) 

11. 

(P3i21,P3221) 

4. 

(P4i2i2,P432i2) 

8. 

(P6222,P6422) 




TABLE I. The 22 chiral space groups in class II, divided in 
11 enantiomorphic space group pairs. Mirroring a crystal be¬ 
longing to one of these space groups will change the crystal’s 
space group to the other one in the enantiomorphic pair. 


chirality of the material, but not its projection on this 
plane. It is clear that in determining the chirality, one 
has to find a way to include the third dimension. The 
most intuitive way of doing this is by rotating the crys¬ 
tal and taking several zeroth order Laue zone (ZOLZ) 
diffraction patterns, as proposed by Goodman and Sec- 
omb [1]. However, in order to determine the absolute 
rotation direction with respect to the crystal from the 
diffraction pattern, one has to perform dynamical sim¬ 
ulations, which require knowledge about the thickness 
of the sample and the structure of the crystal. Taking 
several diffraction patterns under different angles can be 
avoided when looking at the so-called higher order Laue 
zones (HOLZs) in the diffraction pattern [2]. These spots 
contain information about the crystal in the direction 
perpendicular to the projection plane and are suited to 
determine the chirality. However, in the kinematical ap¬ 
proximation Eriedel’s law imposes extra symmetries on 
the Fourier transform of the crystal potential that pre¬ 
vent the chirality from being determined from only one 
diffraction pattern [3]. Therefore, dynamical scattering 
is needed where, again, simulations have to be compared 
with the experiment. An additional difficulty is that not 
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all HOLZ diffraction spots are sensitive to the chirality 
of the crystal and one has to identify these in order use 
them to determine the chirality [4, 5]. 

In previous work [6] we proposed the use of electron 
vortex beams to determine the chirality of crystals with 
a threefold screw-axis. As described theoretically by Nye 
and Berry [7], these electrons’ wavefunction has the form 

= ( 1 ) 


these components in a transmission electron microscope 
(TEM) and discuss the main advantages and obstacles 
of the technique. We start with analytical calculations 
in the first order Born approximation, extend this to in¬ 
clude multiple scattering, and verify our findings with 
multislice simulations. 

II. THEORETICAL FORMULATION 


with r, (f>, and 2 the cylindrical coordinates where the 
z axis is the propagation axis. A vortex electron is 
an eigenstate of the angular momentum operator = 
and possesses a well-defined amount of angular 
momentum of mh per electron, where the number m 
is called the topological charge of the vortex [8]. The 
first experimental realization of optical vortex beams [9] 
inspired active research which led to numerous applica¬ 
tions ranging from nano-manipulation [10-12] and astro¬ 
physics [13-16] to telecommunication [17-19]. Also, elec¬ 
tron vortex beams can be created in a transmission elec¬ 
tron microscope using a phase plate of stacked graphite 
layers [20], holographic reconstruction [21], lens aberra¬ 
tions [22], or the Aharonov-Bohm effect at the tip of a 
magnetized needle [23] , and recent research suggests pos¬ 
sible applications in electron energy-loss magnetic chi¬ 
ral dichroism experiments [21], spin-polarization devices 
[24] , magnetic plasmons [25] and nano-manipulation [26] . 
Since vortex beams are chiral, it may be expected that 
they interact differently with enantiomorphic crystals. 
Indeed, we showed that when focusing a vortex over an 
odd screw-axis with only one heavy atom in its vicin¬ 
ity, the symmetry of the HOLZs will depend on the chi¬ 
rality of the screw-axis [6]. This gives a quantitative 
method to determine the chirality of crystals in trans¬ 
mission electron microscopy, where multiple scattering is 
not required and dynamical simulations can be avoided. 
It also demonstrates the use of electron vortex beams and 
their chiral character in crystallography and shows that 
modifying the electron probe can be of use to investigate 
certain properties. The technique, however, is impeded 
by the extreme sensitivity of the probe position which 
has to be maintained exactly on top of the screw-axis 
during the time of acquisition. Additionally, the symme¬ 
try argument can only be used to determine the chirality 
of odd screw axes. For even screw axes, the diffraction 
pattern still has to be compared with numerical simula¬ 
tions. 

In this work we propose an alternative way of looking at 
electron diffraction. Instead of decomposing the scatter¬ 
ing amplitude in Fourier space as a set of plane waves, we 
write the scattering amplitude as a sum of electron vor¬ 
tices described Eq. (1). Using this framework, we study 
electron scattering on chiral crystals with screw-axis sym¬ 
metry and find that by looking at the 0AM coefficients 
of the scattered wave the chirality of such crystals can 
be determined without the need for dynamical simula¬ 
tions, knowledge about the thickness of the sample, or 
the crystal structure. We also propose a way to measure 


The scattering amplitude for an incident plane wave 
electron with wave vector k to be scattered by a potential 
V{r) to a plane wave with wave vector k' is called the 
scattering amplitude A(fe') = (j)k', K) ^ith 

and k'^ the radial, angular and z component in cylindrical 
coordinates. Since this function is periodic in the polar 
coordinate (pk', it can be expanded into a set of discrete 
Fourier coefficients [27]: 


A{k'j_,(l)k',K) = 








( 2 ) 


In doing so, we write the scattering amplitude as a sum 
of vortices with a well-defined 0AM of mh [see Eq. (1)]. 
The 0AM coefficients am{k_\_, k^) are given by 

1 

am{k^,k,) = — (3) 

V2nJ 0 


and play an important role in the description of the angu¬ 
lar momentum of the total scattering amplitude [28, 29]. 
In Appendix A we show how in the first-order Born ap¬ 
proximation the scattering amplitude for an incoming 
plane wave electron with wave function 4>(r) = to 
scatter to a plane wave 'I''(r) = ’’ on a potential V (r) 

is given by 


A{k') = Y,—Vm{k'^X-k^)e 

— V ZTT 
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with 
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the cylindrical Fourier coefficients of the potential, anal¬ 
ogous to the polar Fourier coefficients in two dimensions 
[27]. Note that these coefficients depend on the origin of 
r, which defines the point around which the expansion is 
taken. Shifting the potential with respect to the origin, 
will result in different coefficients. 

For a crystal with periodicity in the z direction, they can 
be rewritten as 

-I pP poo p2'7T 

Vm,n{k±) dr / d(j)rV{r) 

^ J 0 Jo Jo 

X (6) 
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where P is the periodicity in the z direction and = 
n^, n e N. Because of the periodicity, the latter is 
quantized, which results in a quantized transferred for¬ 
ward momentum kg — k'^ in Eq. (4). 

In what follows we consider consider only elastically scat¬ 
tered electrons for which the energy and the magnitude 
of the wavenumber are conserved during the scattering 
process. We thus have k^ = k'^ -\- k'^ and can write the 
scattering amplitude as 


m 



(7) 


with n = ^{k'^—ko) and k'j_ = i/fep — k'^. Because of the 
discretized transverse momentum, the diffraction pattern 
will consist of discrete rings that can be labeled with the 
number n and are known in conventional electron diffrac¬ 
tion theory as the zeroth order Laue zone (n = 0), the 
first order Laue zone (n = 1), and the higher-order Laue 
zone (n > 0), hereafter referred to as ZOLZ, FOLZ, and 
HOLZ, respectively. Note that in this assumption, the 
ZOLZ would only consist of one single point, while in 
reality it consists of a region of Bragg spots. The reason 
for this is that the crystal has a finite size in the z di¬ 
rection, blurring out the reciprocal lattice points of the 
potential in Fourier space. This effect is mainly visible 
in the ZOLZs where the Ewald sphere lies relatively flat 
with respect to the zone axis, but it is less important for 
the HOLZs where the Ewald sphere cuts the harmonics 
of the z periodicity under an angle [30]. 

When a crystal possesses screw-axis symmetry, it is in¬ 
variant under the combined transformation of rotating 
it around a point on the screw-axis over a certain angle 
A(j) while translating it over a distance Az parallel to the 
rotational axis. Because of the periodicity of the crystal, 
there are only a limited amount of possibilities for 
and Az, 


A(j) = — and Az = —, ( 8 ) 

with M G {2,3,4, 6 } the order of the rotation, P the 
period of the crystal along the screw-axis direction, and 
I G jo,..., M — 1} the translation along the axis in units 
of P/M. In the Hermann-Mauguin notation of space 
groups used in Table I, a screw-axis is denoted as M;. 

In Appendix B we show that the screw-axis symmetry 
puts restrictions on the polar Fourier coefficients in Eq. 
( 6 ) when the crystal is expanded around a point on the 
screw-axis. For a crystal with Mi screw-axis symmetry, 
the nonzero coefficients in the expansion have to suffice, 

m = NM — nl (9) 

with N G'H. From Eq. (7), it is clear that every n-th 
diffraction ring can be seen as a sum of 0AM eigenstates 
^ ^jjg components are given by the coefficient 

Vm.n- As a result, the restriction in Eq. (9) will also be 
reflected in the 0AM states contributing to the different 


Laue zones. 

As an example, consider a crystal that belongs to space 
group P 62 for which M = 6 and 1 = 2. In this case 
the first order Laue zone, n = 1, will only have 0AM 
coefficients for which m G {..., — 8 , —2,4, 10 ,... }, while 
its enantiomorph, belonging to space group P 64 will only 
have components m G {..., —10, —4, 2, 8 ,... }. Here the 
periodicity at which the 0AM components occur is given 
by M. 

It is clear that if we would be able to measure the 0AM 
components of a HOLZ, we could directly identify the 
crystal’s chirality. Moreover, it is sufficient to determine 
for a limited amount of coefficients whether they are zero 
or not, which is generally easier to detect than to mea¬ 
sure their actual value. Note that this property arises 
because of the symmetry and is independent of the crys¬ 
tal structure itself. 


III. MEASURING OAM COEFFICIENTS OF 
HOLZS NUMERICALLY 

In order to investigate the OAM decomposition of scat¬ 
tered electrons in more detail, we perform multi-slice sim¬ 
ulations with the program STEMSIM [31]. We study 
a-quartz that belongs to space group P3i21 or P 3221 , 
where we adopt the structure from Baur [32]. The two 
enantiomorphs posses a right- or a left-handed screw-axis 
for which M = 3, and I = 1 ans I = 2, respectively. In 
Fig. (la) the simulated image is shown for a sample of 
20 nm thickness. 

In Fig. (lb) the diffraction pattern is calculated by tak¬ 
ing the square of the Fourier transform of the simulated 
exitwave. For the ZOLZ, the disk of spots in the middle, 
contains information on the projection of the potential on 
the plane of view. The circle of spots around the ZOLZ 
is the FOLZ and contains information of the crystal po¬ 
tential in the direction perpendicular to the plane of view 
and therefore can be used to determine the chirality of 
the crystal. Since the FOLZ involves high angle scatter¬ 
ing, its signal is considerably weaker than the HOLZ due 
to the atomic scattering factors. 

From the Fourier transform of the simulated exit wave, 
we can numerically calculate the coefficients am,n for any 
n-th Laue zone. First all pixels in Fig. 1(b) are set to 
zero, except for the ones in a small region around the 
Laue zone of interest. Then each pixel is multiplied by 
a factor exp(—imt/)). The coefficients then are obtained 
by adding up all pixel values which is the numerical ap¬ 
proximation of the integral in (3). 

The magnitude squared of the coefficients that build up 
the FOLZ, |am,i(fc_L)P, rn G {—10,10} for right- and 
left-handed a-quartz, sample thickness 20 nm, are shown 
in Fig. (2). Consistent with Eq. (9), it can be seen 
that only the coefficients am,n for which m = NM — nl 
appear in the OAM decomposition. In our example 
for the left-handed a-quartz {I = 2, M = 3) and se¬ 
lecting the FOLZ (n = 1) only the coefficients m = 
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FIG. 1. Simulated TEM image (Left) and diffraction pat¬ 
tern (Right) of right-handed a-quartz with a SOOkeV electron 
beam, samplethickness 20 nm. The circle of Bragg spots in 
the middle of the diffraction pattern is the zeroth order Laue- 
zone that contains information about the projection of the 
potential on the image plane. The ring of spots around the 
ZOLZ is the first order Laue-zone. These spots contain infor¬ 
mation of all three dimensions and thus the chirality of the 
crystal. 



FIG. 2. Simulated GAM coefficients, \am,i{k±)\'^, m G 
{ — 10,10} of the FOLZ in the diffraction pattern of left- (red) 
and right- (blue) handed a-quartz, sample thickness 20 nm, 
assuming the optical axis to be centered on top of a screw-axis. 
Consistent with Eq. (9), for the left-handed enantiomorph 
only the coefficients with m G {... , —2,1,4,...} are present, 
while for the right-handed variant m G {... , —4, —1, 2,... }. 


N.3 — 1.2 G {..., —5, —2,1,4,...} contribute, while for 
the right-handed enanthiomorph the non-zero coefficients 
are the ones for which m G {..., —4, —1, 2,5...}. 

IV. MEASURING OAM COEFFICIENTS OF 
THE HOLZS IN A TEM 

So far we have looked only at OAM decompositions 
by means of simulations where we knew both the 
amplitude and the phase of the scattered wave. In a 
transmission electron microscope, however, we can not 
detect the phase of the scattered wave directly, but we 
can manipulate the scattered wave such that we can 
detect the OAM coefficients, as we demonstrate below. 
Let us consider the image created by selecting a HOLZ 
with an annular aperture in the diffraction plane while 




Projector lens 


FIG. 3. Schematic illustration of experimental setup to mea¬ 
sure the OAM coefficients flAm,!- An incoming plane wave is 
scattered on a sample, after which an annular aperture in the 
objective plane selects only the FOLZ electrons. At the same 
time a spiral phase plate adds an OAM of Am to the FOLZ. 
Using a projector lens, we form an image with the resulting 
electrons where the intensity corresponds to |oAm,iP when 
the decomposition is made around the corresponding point in 
the sample. 


inserting a spiral phase plate as proposed in [33]. This 
experimental setup is schematically shown in Fig. 3. 
With a spiral phase plate we mean any phase plate 
that adds an angular-dependent phase of the form 
to the wave in diffraction space. This can be 
achieved by passing the wave through a material with 
continuously increasing thickness proportional to the 
angular coordinate (j) [20, 34], or by exposing the wave 
to a magnetic monopole-like field induced by the tip of 
a magnetized needle [23]. These spiral phase plates are 
not to be confused with the spiral aperture described in 
[35], which only generates a set of vortex beams with 
different focal points and does not add net OAM to the 
wave. 


Mathematically, the resulting wave in the image plane 
is given by the two dimensional Fourier transform of the 
scattered wave multiplied by a d function selecting the 












































5 


HOLZ and an 0AM eigenstate e 

/ OO n27T 

dk± / d(l)kk±A{k±,(j)k,k^)e~"^'^'^’‘ 

0 Jo 

X S{k_L - (10) 

The point now is that the center of this image (r = 0) 
gives us, up to some constants, the expression for the 
0AM coefficients of the HOLZ [Eq. (3)], 

^Am,n OC 'i/^Am,n(0); 

where again the number n = ^{k'^ — ko) and 
k'^ = ^kl - kf. This means that the intensity at 
the center of the image is proportional to |aAm,n|^- 
Inserting a phase plate imposing an 0AM of —Am to 
the wave can reveal whether or not the Am-th coefficient 
is present in the decomposition around a certain point 
in the crystal simply by looking at the intensity of the 
corresponding point in the image. This is demonstrated 
in Fig. (4), where a screw-axis lies in the center of the 
image. For the left-handed a-quartz, a bright spot can 
be seen on the screw-axis for Am = 1 and Am = 4, 
while for right-handed a-quartz this is the case for 
Am = 2, which is consistent with Eq. (3) and with Fig. 
( 2 ). 

An important note here is that in the theoretical de¬ 
scription of the reciprocal space, we have always put the 
center of our coordinate system on top of a screw-axis. 
When the crystal is shifted, the screw-axis symmetry 
with respect to the center of the coordinate system will 
be broken and the components in the 0AM spectrum no 
longer have to satisfy (9). However, a shift of the sample 
will only result in a shift of the images. The chirality of 
the crystal can therefore still be seen from the intensity 
at the location of the screw axes, which can easily be 
recognized from the three-fold rotation symmetry that 
is present at these points. This is further demonstrated 
in Appendix C for crystals with space groups P4i vs. 
P 43 and P 6222 vs. 6422 . In the ^ 62 ( 4 ) case a clear 
distinction can be made between the sixfold screw-axis 
on the one hand and the threefold screw axes that are 
also present, on the other. This shows that, in principle, 
we can detect all points of screw-axis symmetry in the 
field of view on the atomic level, making it the ideal 
method to study screw-axis symmetry behavior on grain 
boundaries or screw dislocations. 

Because of dynamical scattering, the diffraction pattern 
as well as the 0AM coefficients depend heavily on the 
thickness of the sample and the exact structure of the 
crystal. However, in Appendix D we show that the 
restriction on the components only depends on the 
screw-axis symmetry, independent of the thickness or 
the structure. This means that, in contrast to previously 
developed techniques, when determining the chirality 
from the 0AM coefficients of the HOLZs, we no longer 
require dynamical simulations nor any knowledge about 


left right 



FIG. 4. Comparison of images of left- and right-handed 20 nm 
thick a-quartz. The upper images show the simulated TEM 
image on which no difference can be seen between the enan- 
tiomorphs. The images below are the ones obtained by select¬ 
ing the FOLZ in Fourier space while changing its OAM by an 
amount (in units of h) Am € {0,1, 2, 3, 4, } of the region in¬ 
dicated by the red square. The upper scale bar applies to the 
two upper images, while the lower scale bar applies for the 
OAM -Hltered images. Clear differences can be seen between 
the enantiomorphs of which the appearance of bright spots on 
the location of the screw-axis can be linked to the handedness 
of the screw-axis. 


the crystal structure nor the thickness of the sample, as 
was the case in our previous work using a vortex probe 
to determine the chirality of a threefold screw-axis [ 6 ]. 
The big advantage now is that we can determine the 
chirality of any screw-axis, not only the three-fold, in 
a qualitative way. Additionally, the major drawback. 
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FIG. 5. Simulated OAM coefficients, \am,i{k±_)\'^, m € 
{ — 10,10}, of the FOLZ in the diffraction pattern of left- (red) 
and right- (blue) handed a quartz, sample thickness 20 nm, in¬ 
cluding phonon scattering. Previously forbidden components 
in Fig. 2 now become visible. 


namely, the exact positioning of the probe on top of a 
nanometer-sized screw-axis, has been overcome. 

Because of the symmetry argument and the qualitative 
nature of the signal, the technique seems very promising. 
However, there are some challenges for experimental 
implementation. The first is that it would be prefer¬ 
ential to have a method to switch the total amount 
Am of OAM added in order to be sensitive to multiple 
OAM coefficients. For now, we only have a magnetic 
needle that allows us to study two coefficients with 
opposite OAM. In theory a flipping needle adding OAM 
Am = ±1 suffices to distinguish between 3i and 82 
screw axes, as well as 4i vs. 43 and 61 vs. 65 screw axes, 
since the difference can be seen by studying the 
components. However, this is not the case for the 62 ( 4 ) 
screw-axis. Additionally, the experiment would be a lot 
more conclusive if one would be able to study a series of 
components giving consistent results. The second and 
major obstacle is the fact the feature we are looking for 
in the images, namely, the bright and darker spots, have 
a size of the order of 0.1 A which requires aberration 
correction to work up to 86 mrad or 5A. (Currently, 25A 
is the limit [36, 37].) The reason for this is that the 
images are made out of electrons with high transverse 
momentum and thus high transverse frequencies. Be¬ 
sides putting extreme requirements on the resolution, 
the small size of the features makes the signal very 
sensitive to phonon scattering, where the atoms deviate 
from their equilibrium position. To illustrate this we 
performed multislice simulations, including phonon 
scattering, on a right- and a left-handed quartz sample 
using the experimental Debye-Waller factors from [38] 
obtained at 20 K, the result being shown in Fig. 5. We 
clearly see previously forbidden components appearing 
that obscure the signal, an effect that increases with 
growing sample thickness. For thin samples, however, 
the largest components are still the ones satisfying Eq. 
(9). 


V. CONCLUSION 

By describing electron scattering on chiral crystals 
with screw-axis symmetry in cylindrical coordinates we 
were able to write the scattering amplitude as a sum of 
eigenstates of the orbital angular momentum operator. 
Using the first-order Born approximation, we found that 
the screw-axis symmetry puts restrictions on the OAM 
states building up the HOLZs and that these restrictions 
are dependent on the order and the handedness of the 
screw-axis. We verified this with numerical multislice 
simulations and found this to remain valid even when 
multiple scattering is dominant. This means that the 
screw-axis symmetry and the chirality of the crystal can 
be determined simply by looking at which OAM coeffi¬ 
cients are zero and which are not. 

We propose a method to measure these coefficients by 
making an image after selecting the FOLZ with an annu¬ 
lar aperture in the objective plane, while modifying its 
OAM by an amount Am using a magnetic needle as de¬ 
scribed by Beche et al. [23] . We calculate that the OAM 
components of the FOLZ then are proportional to the 
intensity of the rotation symmetric points in the image 
and show with numerical simulations that the chirality 
can be seen from bright and dark spots in high-symmetry 
points in the image. 

Although the OAM components depend on the thickness 
of the sample and the exact structure of the crystal, we 
show theoretically and numerically that the restriction on 
the components depends only on the screw-axis symme¬ 
try, independent of the thickness or the crystal’s struc¬ 
ture. This means that our proposed method does not 
require dynamical simulations, nor any knowledge about 
the crystal structure or the thickness of the sample. In 
previous work we demonstrated that this was possible for 
crystals with a threefold screw-axis using electron vortex 
beams. The new method, however, extends this for all 
possible screw axes, while the main disadvantages of the 
vortex technique, namely, the exact positioning of the 
probe on the screw-axis, is overcome. On top of that, all 
screw axes in the field of view can be determined at once 
on the atomic scale, while in the vortex setup we could 
only determine the chirality of one screw-axis at a time. 
The main challenge for experimental implementation is 
the required resolution of the order of 0.1 A for a 300- 
keV probe. Although this lies well within the diffraction 
wavelength, it lies beyond the resolution of current state- 
of-the-art microscopes. Reaching these may require tech¬ 
nological breakthroughs which pose serious challenges. 
Also, adding a range of different amounts of OAM to the 
scattered wave in Fourier space such that we can acquire 
a series of images with a different OAM filter would be 
preferential, where currently we can only impose only two 
values with opposite sign using a single magnetic needle. 
Finally, this work shows that looking at the OAM decom¬ 
position of a scattered wave can reveal extra information 
about a sample if one understands how the scattering 
process influences the OAM of the wave. The experi- 









7 


mental method to measure the 0AM components pro¬ 
posed here would also work without FOLZ filtering, in 
which case it would be directly applicable in a conven¬ 
tional TEM without the extreme performance which is 
needed here. Such a method could find use way beyond 
the field of chiral space group determination. 
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Appendix A: Diffraction in cylindrical coordinates 

Consider a plane-wave electron with wave vector k 
scattered on a potential V{r). In the first-order Born 
approximation, the scattering amplitude to scatter to a 
plane wave with wave vector k' is given by [30] 

A[k')={k'\V\k) 

=Jdre-^’^'-^Vir)e^^r (Al) 

In what follows it will be convenient to have an expres¬ 
sion for the scattering amplitude in cylindrical coordi¬ 
nates. To calculate this, we first expand the potential in 
its cylindrical Fourier components [27], 

/+O0 pOO OQ 

dk^ / dA:_Lfc_L ^ Vm{k±,k:,) 

“'0 m=-oo 

X (A2) 

V ZTT 


Using the Jacobi-Anger identity. 


cos(0-0'j,) 


= (A5) 
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we get 
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dr Jm{k±r)Jm'(.k'^r) 
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(A6) 

The integrals over r, z, and (j) give two Dirac 5 functions 
and a Kronecker 6. The final expression for the scattering 
amplitude thus becomes 


A(fc') = 


-| ^ + oo poo ^ 

') =^= / dkz / dfc_LA:_L Yl K^(fc_L,A:z) 

V 27r J —oo Jo / 

* m m' = — f-v~i 


m.m — — OO 


X I 
1 


+ ko- k',) 

fe_L 

i-^Vrn{k'^X - ko)e ^^‘^^'. (A7) 


Equation (A7) gives the most general form of a scat¬ 
tering amplitude for a plane wave electron to scatter to 
a plane wave with wave vector k' scattered on a po¬ 
tential determined by its cylindrical Fourier coefficients 
Vm{k'j_,k'^ - ko). 


with 


/ +00 poo p 2 n 

dz dr d(j)rV{r) 

— OO Jo JO 


X (A3) 
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For the scattering amplitude for an incoming electron 
along the z axis with wave number k^, substituting Eq. 
(A2) into (Al) gives us 


Appendix B: Scattering on potential with screw-axis 
symmetry 

For a crystal potential with screw-axis symmetry 
around the origin one has 

V{r,(j},z) = V{r,cj) + A(l),z + Az), (Bl) 

where the periodicity of the crystal puts restrictions on 
A(j) and Az, 


/ ■ /*+oo poo oo 

dr dkj dk^kj_ ym{k±,k,) 

J Jo m=-oo 






(A4) 


A(j) = — and Az = —, (B2) 

M M ^ ' 

with P the period of the crystal along the screw-axis 
direction, M G {2,3,4, 6} the order of the rotation and 
I G (0,..., Af — 1} the translation along the axis in units 
of P/M. 










A first thing to note is that because of the periodicity in 
the z direction, the potential can be written as 


/ OO 

dk^ fcx V Vm,n{k±) 
0 _ _ 


m,n— — oo 


X (B3) 
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with 


and the relation fcg = + k'^ has to be satisHed. Since 

for a given k'^ or equivalently a number n, the variable 
k± is fixed by the relation k± = ^Jkl — k'^ and can be 
dropped in Eq. B 8 . We then get 


-dn ) 
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(B9) 


where 
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(B4) 


Secondly, the screw-axis symmetry puts restrictions on 
the contributing components. 

2 oP poo ^27r 
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^ J 0 Jo Jo 
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Vjn,n - (BIO) 

The discretization of the transverse momentum makes 
that the diffraction pattern consist out of discrete rings 
which can be labeled with the number n. These rings co¬ 
incide with the ZOLZ (n = 0), FOLZ (n = 1) and higher 
order Laue zones (n > 0) known from conventional elec¬ 
tron beam diffraction theory. 

The scattering amplitude for an incident plane wave elec¬ 
tron with wave vector k to be scattered by a potential 
V{r) to a plane wave with wave vector k' is called the 
scattering amplitude A{k') = A{k'j_, , k'^) with k'j_,(j)i^' 

and A:(, the radial, angular and z component in cylindrical 
coordinates. Since this function is periodic in the polar 
coordinate (j)k>, it can be expended into a set of discrete 
Fourier coefficients [27]. 

We can write the scattering amplitude in cylindrical 
coordinates A{k') = A{k'j_,(j)k',k'^) with k'j_,(j)k' and k’^ 
the radial, angular and z component. Since this function 
is periodic in the polar coordinate (pk', it can be expended 
into a set of discrete Fourier coefficients [27]. 


A{r^,cPk',ki) = Y.am(.k'^,k',)e^^t (Bll) 

V ZTT 


Therefore the non-zero components have to satisfy the 
relation 

Stt 

mA(j) + n — Az =N2Tr (B 6 ) 

or equivalently 

m = NM — nl (B7) 

with N G N. Polar Fourier coefficients for which this 
relations is not fulfilled, will automatically be zero. 
Knowing this, we can write the scattering amplitude, Eq. 
(A7), for a potential with screw-axis symmetry as 


In doing so, we write the scattering amplitude as a sum 
of vortices with 0AM mh, see Eq. (1). The 0AM coef¬ 
ficients am{k±, kz) are given by 

1 7^’^ 

arn{k^,kz) = ^ dcPkA{k^,Pk,kz)e-^"^^K (B12) 

The use of cylindrical coordinates makes it clear how ex¬ 
actly the scattering amplitude is written in terms of 0AM 
eigenstates. From Eq. (A7) it follows that these are pro¬ 
portional to the cylindrical expansion coefficients 

= -*(,). (B13) 


A(fc') 








E --NM^r (p \J{NM-nl)(j>y 

^ yNM-nl,n(Kj_)e 

N 


(B 8 ) 


For a crystal with screw-axis symmetry we thus get, Eq. 
(B7) and (B9) 


^m,n 


i 'm = NM — nl 

0 otherwise 


(B14) 


were the transferred forward momentum, k'^ — ko = 
is discretized because of the periodicity of the potential. 
For elastically scattered electrons the energy is conserved 


As an example consider a crystal that belongs to space 
group P 62 , for which M = 6 and 1 = 2. In this 
case the 0AM decomposition of the first order Laue 
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zone, n = 1, will only have components for which 
TO G — 8 , — 2 ,4, 10 ,... }, while its enantiomorph, 

belonging to space group P 64 will only have compo¬ 
nents TO S — 10 , —4, 2 , 8 ,... }. M allways gives 

the periodicity of the occuring components, while —I, 
in the FOLZ-case, is given by the largest negative 
TO-component contributing. 


Appendix C: Examples of crystals with other 
screw-axis symmetries 

So far we only looked at a quartz with a threefold 
screw-axis. Of course relation (B14) puts restrictions on 
the 0AM components of HOLZs of diffraction patterns 
of crystals belonging to all chiral space groups in Table 
I. This is demonstrated in Fig. ( 6 ) where we look at 
the 0AM decomposition of P 6222 and P 6422 (3 quartz 
[39] with respect to the sixfold screw-axis and P4i and 
P 43 a-tricesium heptaphosphide [40]. It is clear that in 
all cases, indeed, only the 0AM components satisfying 
relation (B14) are present in the FOLZ. 




FIG. 6 . Simulated OAM coefficients, \am,i{k±3)\^, rn € 
{ — 10,10}, of the FOLZ in the diffraction pattern of left- (red) 
and right- (blue) handed (a) CS 3 P 7 and (b) j3 quartz.Sample 
thickness is 20 nm and 200 nm, respectively. Again, relation 
(B14) is fulfilled. 

A big advantage of our technique is that one can iden¬ 
tify different screw axes on one series of images, as can 
be seen in Fig. (7). Here we look at the FOLZ-filtered 
image of P 6222 /3-quartz with different Am. As expected 


from Fig. ( 6 b), we can identify a 62 screw-axis from the 
bright spots at Am G {4,10}. However, at the same 
time we can see bright spots in the threefold symmetric 
points at Am = {1,4,7,10,12} indicating a 3i screw- 
axis, which also is present in space group P 6222 . The 
technique is thus capable of determining all screw axes 
and their handedness on the atomic scale in the entire 
field of view. 

In Fig. 8 we give the FOLZ-filtered images for tricesium 
heptaphosphide. Because of the relatively large period¬ 
icity in the z direction, the FOLZ will have a smaller 
radius. This means that the electrons that make up the 
image will have a smaller transverse frequency, which is 
why the features now are of the order 0.5 A wide, a factor 
of 5 bigger than previous examples. Lowering the accel¬ 
eration voltage will also result in a smaller radius for the 
FOLZ, having the same effect on the size of the spots we 
are trying to distinguish. However, doing this generally 
results in a lower resolution of the microscope. 



FIG. 7. TEM image and FOLZ-filtered, OAM modified im¬ 
ages of right-handed /3 quartz, Am € {1,4,7,10,13} . From 
the bright spots a 62 screw-axis in the center and six 3i screw 
axes at the edges can be identified. Scale bar applies to all 
subfigures. 


Appendix D: Effect of multiple scattering 

Of course, electron scattering is mostly dynamical, 
making the FOLZ and its OAM spectrum dependent on 
the thickness of the crystal. The latter can be seen in 
Fig. (9), where the OAM spectrum is shown for left- 
and right-handed a quartz for thicknesses ranging from 
50 nm to 150 nm. However, from this hgure it appears 
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FIG. 8. TEM image and FOLZ-filtered, OAM modified im¬ 
ages of right-handed CS3P7, Am € { — 1,3, 7} . From this a 4i 
screw-axis in the center and in the corners can be identified. 
As a result of the relatively large periodicity in the direction 
of the screw-axis, the FOLZ will have a smaller radius, which 
results in bigger features in the image of the order of 0.5 A. 
The scale bar applies to all subfigures. 


that all forbidden coefficients in Eq. (B14) remain 
forbidden even when multiple scattering is dominant, 
independent of the thickness of the sample. 

In order to understand this, we first look at the effect 
of a translation in the z direction and a rotation around 
the optical axis of the incoming beam on the scattering 
amplitude, as illustrated in Fig. (10). The effect of a 
translation of the sample along the z direction, the prop¬ 
agation direction of the incoming electron, can be fully 
incorporated simply by adding a phase to the incom¬ 
ing electron. As an effect of this, the scattering amplitude 
of a shifted potential V'{r,(j),z) = V{r,4>,z + Az), also 
gets this extra phase such that the scattering amplitude 
is written as 

fc.) = (Dl) 

Since the incoming electron wave is rotationally invari¬ 
ant, rotating the sample simple results in a rotated scat¬ 
tering amplitude. This means that, when rotating the 
sample along the optical axis over an angle A(j>, the scat¬ 
tering amplitude now becomes 

A'{k±,(j)k, k^) = A{k±, (j)k + A(j), fcz). (D2) 

When combining a rotation of the sample with a trans¬ 
lation in the z direction, V'{r) = V{r,(j) + Acj), z + Az), 
we thus get for the scattering amplitude 


A'ik^, cj>k, k,) = Aik^,<j,k + z)e*'=^^A (D3) 

It is important to note that the treatment includes full 
dynamical scattering as no kinematical approximation is 
made. 




FIG. 9. Simulated OAM coefficients, |a„,,i(fcx)p, m G 
{ — 10,10}, of the FOLZ in the diffraction pattern of left- (red) 
and right- (blue) handed a quartz for sample thicknesses (a) 
50, (b) 100 and (c) 150 nm. Although the coefficients depend 
on the thickness of the sample, they will still obey relation 
(B14), even when multiple scattering is dominant. 


Let us now consider the special case of a crystal with 
screw-axis symmetry for which the potential has the 
property 


V{r,(j), z) = V{r,(j) + A(j), z + Az) (D4) 

for A(j) and Az, as before 

A(p = — and Az (D5) 

M M ^ ' 

with M S {2,3,4, 6} the order of the rotation, P the 
period of the crystal along the screw-axis direction, and 
I G (0 ,..., M — 1} the translation along the z axis in 
units of P/M. 

Since the crystal potential is invariant under this trans¬ 
formation, the scattering amplitude has to be so as 
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Also, in the dynamical approximation, conservation of 
energy is assumed, such that we can rewrite the scat¬ 
tering amplitude as A„((/)) with n = ^{k'^ — ko), k'j_ = 
~ ^'z- this notation the Eq. (D6) becomes 


^n{4>k) — An ( (t>k + ^ ) 6* " ^ (D7) 


As before, we can calculate the 0AM expansion coeffi¬ 
cients for the scattering amplitude using the formula 


1 

—d(/)feAn((Afc)e-™‘^^ (D8) 



Filling in relation (D6), eventually gives 


^m,n 


1 /*27r 

= ^y ^ d^k' An 


(D9) 

(DIO) 


Therefore the non-zero components have to satisfy the 
same relation as before 


FIG. 10. Sketch of the effect of translating a sample along 
the beam direction and a rotation of the sample around the 
optical axis. The first adds an extra phase to the scat¬ 
tering, while the second results in a rotation of the scattering 
amplitude. 

well. This means that in Eq. (D3) A'(k±,(j)k,kz) = 
A{k±, (j>k, kz) and thus 

A{kj_,(j)k, kz) = A{k±,(j)k + Acjife, (D6) 


mA(j) + n — Az =N2 tt, (DU) 

or equivalently, 

m = NM-nl. (D12) 

Therefore the nonzero 0AM components of the scatter¬ 
ing amplitude have to satisfy the relation (9), and in the 
case of dynamical scattering also. 
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